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PML  1994-A33 
A92/KM/407 


In  het  kader  van  de  opdracht  A92/KM/407  “Kwetsbaarheid  van  anti-schip  geleide  wapens”  is 
een  code  ontwikkeld  welke  een  analyse  beoogt  te  kunnen  geven  van  het  gedrag  van  een  inko- 
mend  geleid  wapen  dat  wordt  bestreden  door  verdedigingsmiddelen  van  het  schip  (of  vloot)  dat 
wordt  aangevallen.  Met  deze  code  kunnen  de  schadelijke  gevolgen  van  een  exploderende 
warhead  op  een  inkomend  anti-schip  geleid  wapen  worden  bepaald.  In  het  rapport  wordt  een 
algemene  beschrijving  gegeven  van  de  Missile  ^sterns  Damage  Assessment  Code  (MISDAC), 
die  ontwikkeld  werd  uit  een  tweetal  bestaande  codes  namelijk  Target  Vulnerability  Assessment 
Code  (TARVAC)  en  Missile  Vulnerability  Assessment  Code  (MISVAC). 

De  code  berekent  voor  iedere  onderschepping  welke  onderdelen  van  het  geleide  wapen  worden 
beschadigd,  waarbij  rekening  wordt  gehouden  met  de  onderscheppingsgeometrie,  de  relatieve 
snelheden,  de  verschervende  werking  van  de  warhead,  en  de  penetratie  van  de  scherven  in  het 
doel.  De  locatie  van  de  springpunten  in  de  buurt  van  het  doel  kan  enerzijds  worden  bepaald  met 
behulp  van  speciale  programmatuur,  of  anderzijds  via  geleverde  springpunten.  Ook  is  het  moge- 
lijk  de  schade  door  projectielen  te  bepalen,  bijvoorbeeld  afkomstig  van  snelvuur  kanonnen 
(Goalkeeper). 

De  reden  voor  de  ontwikkeling  van  MISDAC  is  het  verkrijgen  van  een  beter  inzicht  in  de  be- 
schadigingen  die  kunnen  ontstaan  in  het  doel,  eventueel  gespecificeerd  naar  type  onderdeel,  om- 
dat  in  de  voormalige  code  MISVAC,  die  tijdens  de  NATO  Anti-Air  Warfare  System  (NAAWS) 
periode  werd  ontwikkeld,  slechts  grotere  delen  van  het  doel  (resp.  secties  of  sectoren)  konden 
worden  geanalyseerd.  Hierdoor  is  MISVAC  beperkt  tot  de  drie  in  NAAWS  gedefmieerde  doelen 
vanwege  de  door  de  Britten  ingebrachte  “experimentele”  gegevens.  De  feitelijke  achtergronden 
zijn  onduidelijk  en  maken  van  dit  model  dan  ook  een  “black  box”  benadering. 

Ten  opzichte  van  de  langer  bestaande  code  TARVAC  bestaat  nu  de  mogelijkheid  om  schade  ten 
gevolge  van  directe  treffers  op  afzonderlijke  componenten  nader  te  analyseren,  alsmede  de 
schade  door  blast.  Laatstgenoemde  mogelijkheid  is  echter  op  eenvoudige  wijze  geimplemen- 
teerd,  omdat  een  dynamische  berekening  van  constructies  niet  in  een  schootslijnmethode  kan 
worden  uitgevoerd. 

Nadat  de  schade  berekeningen  zijn  uitgevoerd  wordt  door  een  andere  code  Flight  Path 
Emulation  of  Missiles  (FLIPSIM)  bepaald  of  het  beschadigde  missile  in  staat  is  het  schip  alsnog 
te  treffen,  gegeven  de  interceptie  afstand  (ook  wel  beschadigingsafstand  genoemd)  en  het  type 
manoeuvre  dat  wordt  uitgevoerd.  De  totale  uitkomst  van  de  berekeningen  geeft  de 
uitschakelkans  van  het  missile,  gegeven  een  (succesvolle)  onderschepping.  Tevens  bestaat  de 
mogelijkheid  om  op  basis  van  het  gedrag  van  het  beschadigde  missile  te  komen  tot  bepaalde 
klassen  van  (verstoorde)  banen,  hetgeen  een  hulpmiddel  kan  zijn  voor  de  zogenaamde  “kill 
assessment”  studies. 
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In  het  rapport  wordt  ingegaan  op  vele  details  met  betrekking  tot  de  schootslijn  benadering,  de 
implementatie  van  blast  effecten  en  de  wijze  waarop  de  diverse  uitschakel  mogelijkheden  wor- 
den  verwerkt. 

Het  rapport  concludeert  met  het  feit  dat  MISDAC  is  ontstaan  uit  de  twee  oudere  codes,  namelijk 
de  “standaard”  code  TARVAC  en  de  daaruit  afgeleide  code  MISVAC.  Hierbij  zijn  van  beide  co¬ 
des  de  nadelen  zoveel  mogelijk  vermeden,  en  de  voordelen  benut. 

De  uiteindelijke  kans  op  uitschakeling  van  het  missile  moet  nog  worden  gei'mplementeerd  via 
geavanceerde  zogenaamde  6-DOF  simulaties  met  behulp  van  FLIPSIM.  De  koppeling  van 
FLIPSIM  aan  MISDAC  is  een  aspect  dat  in  het  vervolgonderzoek  plaats  zal  vinden. 
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Amin  ■  minimum  presented  area 

Amax  -  maximum  presented  area 

Ap,  Ap  -  (average)  presented  area 

Aw  -  wetted  surface 

A,  B,  C  Westine  constants 
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Cd  -  drag  coefficient 
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D  -  (projectile)  diameter 
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m  -  projectile/fragment  mass 

M  -  Mach  number:  M=V/ca 
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Pdh  -  probability  of  damage  given  a  hit 

Pkd  -  probability  of  kill  given  damage 
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p  -  pressure 
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R  -  range  between  missile  and  target 

Ra  -  absolute  gas  constant 
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Rc  -  charge  radius  (TNT) 
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r  -  radius 


Sn,  Sn  -  (average)  shape  number 

t  -  time ;  thickness 

tmax  ■  max.  plate  thickness 

tp  -  positive  phase  duration 

At  -  difference  in  time-of-arrival 

T  -  temperature 

u  -  flow  (particle)  velocity 

0,  Ug  -  (average)  shock  wave  velocity 

V  -  velocity 

Vf  -  fragment  velocity 

Vfa  -  average  initial  fragment  ejection  velocity 
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Vfo  - 
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impact  fragment  velocity 

Vfmt  - 
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Vfta  - 

average  relative  fragment-target  velocity 
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Vfti  - 
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Vm  - 

missile  velocity 

Vp  - 

projectile  velocity 

Vpt  - 
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Vt 

target  velocity 

Vr 

residual  velocity 

w 
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W 

warhead  mass,  TNT-eq. 

x,y 

co-ordinates;  variables 

z 

scaled  range  Z=RyW^^^ 

Greek  symbols 
a 

V 

e>ex,y  - 

(p 

P 

% 

Y 

e 


Q 

P>  Pm  ■ 

Pa 

CO 


azimuth  angle 
elevation  angle 
dispersion  angles 

miss  orientation  angle;  angle  of  incidence 

pitch  angle;  Weibull  parameter;  time  factor 

fragment  ejection  angle 

yaw  angle,;  ratio  of  specific  heats,  y=\A 

fragment  normal  impact  angle 
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random  (0,1)  variable 

density  (of  fragment) 

air  density  (sea  level) 

fragment  rotation  speed 
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Obsolete  expression  to  indicate  an  ASM  but  not  to  be  confused  with  Air-to-Surface.  Air-to-Surface  is 
designated  as  Air-to-Ground  (AGM)  for  land-based  missiles  and  as  Anti-Ship  for  naval  versions. 
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1  Introduction 


This  report  has  been  written  to  assist  the  Royal  Netherlands  Navy  in  developing  adequate  means 
to  defend  their  battleships  against  future  threats,  particularly  from  incoming  antiship  missiles.  In 
addition  to  research  being  conducted  on  survivability  enhancement  of  battleships,  the  Directorate 
of  the  Royal  Navy  Materiel  (DMKM)  required  an  investigation  into  the  effectiveness  of  the  air 
defence  systems  currently  installed,  and  to  be  installed  on  board  future  battleships,  within  the 
framework  of  A92/KM/407,  “Vulnerability  Analysis  of  Anti-Surface-Ship-Missiles”.  This  re¬ 
port,  however,  only  deals  with  the  methodology  that  has  been  developed  to  assess  the 
vulnerability  of  missile  targets. 

In  the  Weapon  Effectiveness  research  group  of  the  Weapons  &  Platforms  Division  at  the  TNO- 
Prins  Maurits  Laboratory  (PML),  a  target  vulnerability  assessment  code  (TARVAC)  has  been 
developed.  This  code  assesses  the  effects  of  projectiles  and/or  fragments  with  respect  to  the  ca¬ 
pability  of  killing  a  target.  Basically,  the  method  employs  a  shot  line  technique,  which  represents 
the  trajectory  of  an  actual  projectile  or  fragment.  Several  parts  of  the  code  deal  with  geometrical 
intersections  of  such  shot  lines  with  the  target.  In  the  code,  the  target  is  represented  by  a  com¬ 
puterised  description,  where  components  are  represented  by  solids  (geometrical  figures),  to 
which  certain  material  properties  (e.g.  material  density,  dimensions)  are  assigned. 

When  dealing  with  fast  incoming  antiship  missiles  the  defensive  capabilities  on  board  the  ship 
are  tested  to  the  limits.  The  missile  must  be  engaged  far  enough  away  to  prevent  the  missile  or  its 
debris  from  hitting  the  ship.  If  any  part  of  the  missile  hits  the  ship,  it  may  still  inflict  severe 
damage  on  the  ship.  Apart  from  a  fast  reacting  air  defence  asset,  the  kill  of  the  missile  must  be 
ensured.  To  this  end  all  damaging  effects  having  a  certain  kill  capability  are  addressed  into  more 
detail,  not  only  fragment  impacts,  but  also  direct  hits  (from  a  surface-to-air  missile)  and  blast 
damage  leading  to  structural  kill.  In  this  respect  it  was  found  that  TARVAC  was  inadequate,  and 
hence  an  improvement  was  necessary. 

The  direct  hit  option  is  dealt  with  by  including  so-called  missile  representing  points,  and  re¬ 
garding  relative  missile-target  trajectories  as  equivalent  shot  lines.  The  blast  damage  option, 
however,  cannot  be  dealt  with  using  a  shot  line  approach.  Although  several  options  exist  to 
model  blast  damage,  it  was  decided  to  keep  the  level  of  detail  for  blast  modelling  as  simple  as 
possible.  To  this  end,  so-called  blast-critical  components  were  introduced  in  addition  to  specified 
damage  radius  criteria.  Complex  blast  wave  -  target  interactions  are  then  ruled  out  by  definition, 
but  it  is  emphasised  that  fragment  damage  predominates  the  overall  damage  to  the  target. 
Therefore  a  new  code  has  been  developed  which  in  comparison  with  the  former  code  is  named 
Missile  Systems  Damage  Assessment  Code  (MISDAC).  MISDAC  is  a  further  development  of 
the  Missile  Vulnerability  Assessment  Code  (MISVAC),  which  was  derived  from  TARVAC  and 
the  code  SHAMDAT  (Sectional  Hardness  Methodology  for  Damaged  Air  Targets)  from  the 
UK  [1].  MISVAC  has  been  used  exclusively  in  the  NATO  Anti-Air  Warfare  System  (NAAWS) 
evaluation,  and  has  been  applied  in  the  SEAROADS  project  to  assess  the  lethality  of  the  NATO 
Sea  Sparrow  missile  warhead  [2].  Although  MISVAC  includes  the  options  mentioned  before,  it 
lacks  the  detailed  analysis  of  internal  components  being  damaged,  while  the  overall  kill  assess¬ 
ment  includes  the  target-range-to-go  as  a  simple  factor.  Target-range-to-go  is  defined  as  the 
range  from  the  ship  being  attacked,  at  which  damage  to  the  target  occurs,  also  called  the  intercept 
range.  Although  by  far  one  of  the  fastest  computational  methods,  much  of  the  background  is 
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seemingly  based  on  experimental  data  and/or  educated  guesses  making  the  whole  concept  a 
black-box  like  approach.  MISDAC  is  constructed  in  such  a  way  that  detailed  damage  assessment 
of  components  is  performed  -  similar  to  TARVAC  -  and  that  post-damage  flight  trajectories, 
given  damage  to  the  missile  systems,  are  conducted  by  a  code  called  Flight  Path  Emulation  of 
Missiles  (FLIPSIM)  to  find  the  ultimate  probability  whether  the  missile  or  its  debris  will  hit  the 
ship,  which  is  a  measure  of  the  probability  of  missile  kill.  In  other  words,  the  output  of  MISDAC 
serves  as  input  for  FLIPSIM.  Figure  1  illustrates  the  development  of  MISDAC  showing  the  es¬ 
sential  modules  embedded  in  the  program. 


MISVAC 


MISDAC 


trajectories 


TARVAC 


shotline 

generation 


target 

intersection 

damaged 

components 


killed 

components 

killed 

systems 


Figure  I  Schematic  development  of  MISDAC 


In  the  BURST/GRID  module,  the  shotlines  from  a  detonating  warhead  are  generated  and  the 
fragments  hitting  the  target  perimeter^  are  determined.  In  the  TARIN  module,  a  more  accurate 
determination  is  performed  in  which  components  of  the  target  are  intersected.  In  DAMCO,  all 
components  are  determined  that  are  penetrated  using  ballistic  penetration  relationships.  In 
KILCO,  the  probability  of  kill  of  the  damaged  components  are  determined  using  kill  criteria,  and 
finally  in  KILSY,  the  probability  of  kill  of  all  systems  making  up  the  complete  target  are  deter¬ 
mined  using  special  logical  combinations  according  to  a  kill  tree.  In  MISVAC,  the  KILCO 
module  has  been  replaced  by  the  British  SHAMDAT  module,  while  DAMCO  has  been  stripped 
because  penetrations  are  not  calculated.  An  extra  feature  in  MISVAC  (and  also  MISDAC)  is  the 
possibility  to  use  a  fuse  model  or  a  generic  model  to  create  burst  points.  However,  at  present  the 


The  whole  target  is  enclosed  in  a  box. 
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fuse  model  only  woiks  for  the  three  existing  NAAWS  targets  (i.e.  SBS,  SSS  and  SHD  missiles) 
for  each  of  which  a  so-called  stick-cone  target  model  exists. 

The  computations  in  FLIPSIM  are  performed  using  dynamic  flight  trajectory  simulations,  where 
the  missile  aerodynamics  and  controls  are  incorporated  in  a  six-degrees-of-freedom  (6-DOF) 
simulation  model.  So,  an  assessment  has  to  be  performed  about  the  residual  missile's  control  and 
aerodynamics  properties  after  being  engaged,  which  forms  the  most  difficult  part  of  the  study. 
The  final  outcome  of  the  computations,  i.e.  the  probability  of  kill  of  the  missile  given  an 
engagement  (either  by  guns  or  missiles)  is  then  found  using  a  logical  combination  of  MISDAC 
with  FLIPSIM.  This  report  will  describe  the  basic  features  of  only  the  first  part,  i.e.  MISDAC. 


2  Description  of  the  method 


Basically,  the  method  employed  in  MISDAC  is  quite  similar  to  the  original  method 
TARVAC  [3].  This  method  has  been  developed  over  many  years,  and  changes  have  been  incor¬ 
porated  ever  since.  In  principle,  the  target  is  geometrically  represented  by  seven  solids,  viz.:  box, 
pentahydron  (5-planes),  cylinder,  (truncated)  cone,  parallelepiped,  triangular  plate  (with  thick¬ 
ness)  and  sphere.  To  model  a  component,  several  solids  can  be  used,  in  addition  to  “negative” 
solids,  by  which  voids  etc.  can  be  modelled.  In  the  target  database,  the  components  are  indicated 
by  two  6-digit  numbers  which  have  a  special  meaning.  The  first  6  digits  indicate  the  component 
number,  which  may  be  arbitrarily  chosen.  The  second  6-digit  number  contains  the  component's 
physical  information.  The  first  digit  must  be  either  a  0,  1  or  2,  indicating  its  vulnerability 
(O=non-vulnerable,  1  =vulnerable  to  fragments,  2=vulnerable  to  blast).  The  next  2  digits  represent 
the  material  code  (ranging  from  1  to  25)  and  the  last  3  digits  represent  the  void  ratio  in 
percentage  (varying  from  0  to  100  per  cent). 

One  of  the  differences  between  MISDAC  and  TARVAC  lies  in  the  vulnerability  code  for  a 
component  being  vulnerable  to  blast  (=2).  In  essence,  a  blast-vulnerable  component  is  repre¬ 
sented  by  a  cylinder  with  zero  radius  (or  line  element),  so  that  its  structure  fits  well  in  the  stan¬ 
dard  target  construction  method  using  (the  same)  solids.  However,  shotline  intersections  with  the 
actual  component  are  omitted  since  the  vulnerability  code  (=2)  is  tagged  to  skip  this  kind  of 
operation. 

In  the  missile-target  engagement,  where  an  incoming  sea  skimming  missile  target  is  intercepted 
by  a  surface-to-air  missile  (SAM)  defending  the  ship,  the  engagement  geometry  is  outlined  in 
Figure  2.  The  main  parameters  driving  the  engagement  are: 

•  missile  velocity 

f  which  defines  the  relative  missile-target  velocity  Vmt 

•  target  velocity  Vj  J 

•  azimuth  angle  a 

•  elevation  angle  \(r 

•  miss  distance  MD 

•  miss  distance  orientation  tp 

The  missile  has  a  pitch  (P)  and  yaw  angle  (y)  relative  to  its  flight  path.  Due  to  symmetry  the  roll 
angle  is  ignored. 
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Figure  2  Missile-target  engagement  geometry 

To  maintain  similarity  with  TARVAC,  the  target  modelling  uses  the  same  co-ordinate  system, 
but  the  target  velocity  is  directed  to  the  -x  axis  (left).  The  target  and  SAM  missile  are  represented 
by  their  reference  points  TRP  and  MRP,  which  are  located  exactly  in  the  geometric  centres  of  the 
target  and  missile,  respectively.  Further  explanation  of  the  intercept  geometry  can  be  found  in 
[1].  For  each  engagement,  a  set  of  relative  missile-target  trajectories  is  created  defined  by  a 
(fixed)  miss  distance  MD  and  a  number  of  miss  distance  orientation  angles  (cp).  In  this  way  a 
number  of  parallel  trajectories  are  created  lying  upon  the  cylinder  around  the  target  with  the  miss 
distance  MD  as  its  radius.  Along  each  of  the  trajectories,  a  number  of  detonation  points  are 
selected,  for  example  spaced  at  constant  time  intervals  before  and  after  t=0^.  The  point  t=0  is 
exactly  the  point  where  the  distance  between  missile  and  target  is  equal  to  the  miss  distance.  This 
point  is  usually  referenced  as  Closest  Point  of  Approach  (CPA),  see  Figure  3. 


Figure  3  Shot  lines  -  missile  target  encounter 

a)  -  general  view  b)  -  2-D  view 


1 


Or  one  single  fuse-determined  burst  point. 
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The  results  of  the  computations  are  then  averaged  over  all  burst  points  to  give  an  average  for 
this  miss  distance.  This  generated  encounter  file,  containing  the  possible  detonation  positions  of 
the  warhead  relative  to  the  target,  can  also  be  used  in  other  external  programs  to  estimate  the  ac¬ 
tual  fuse-triggered  detonation  points.  In  fact,  this  has  been  employed  in  earlier  computations  for 
the  Royal  Netherlands  Navy  project  SEAROADS  [2].  However,  data  pertaining  to  warhead  burst 
points  driven  by  proximity  fuses  are  generally  obtained  from  the  TNO  Physics  and  Electronics 
Laboratory. 


3  ProjectileAvarhead  damage  effects 


3.1  Shot  line  generation 

To  analyse  the  effects  of  projectiles  hitting  the  target,  or  the  effects  of  a  warhead  detonation  in 
the  vicinity  of  the  target,  the  shot  line  technique  is  used.  Each  projectile  or  fragment  trajectory  is 
simulated  by  its  corresponding  shot  line,  having  an  identification  number,  a  specified  mass,  ve¬ 
locity,  origin  and  direction.  A  few  options  exist  by  which  a  grid  of  shot  lines  can  be  created.  One 
option  is  to  use  a  parallel  grid  of  equally-spaced  shot  lines,  “fired”  from  a  certain  direction  to  the 
target.  For  statistical  purposes,  a  random  position  within  a  grid  cell  may  be  chosen.  In  the  pro¬ 
jected  plane,  the  shot  line  co-ordinates  are  given  by: 

X-Xn  -l-Q(Xn+,  -x„) 

y~ym  "*’^(ym+i  ~ym) 

where  Q  is  either  fixed  (^2=0.5)  or  a  (0,1)  random  variable  in  an  (n,m)  grid,  see  Figure  4. 


In  another  option,  a  so-called  “aiming”  grid  can  be  selected,  which  is  composed  of  an  aiming 
point  and  a  number  of  shot  lines  “fired”  at  this  aiming  point  using  a  statistical  error,  usually  a 
circular  error  normally  distributed  around  the  aiming  point,  see  Figure  5.  In  the  projected  plane 
the  shot  line  co-ordinates  are  found  using  bivariate  normal  deviates  according  to  Box-Muller  [4]: 
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1/9 

x  =  [-21nQi]  cos(27l^^2) 

y  =  [-2  In  Q|  ^  ^  sin(2jrQ2 ) 


(2) 


where  Q|,  Q2  are  independent  (0,1)  random  variables.  The  (x,y)  deviates  given  a  range  R  and  an 
angular  error  8  yield  the  normal  “hit”-distribution  around  the  “aiming”  point  (0,0): 


1/9 

x^REx  [-21nQ,]  cos(2jiQ2) 
y  =  REy  [-21nQ,]‘'^sin(27iD2) 


(3) 


where  8x  and  Ey  are  two  independent  errors.  If  they  are  equal,  a  circular  error  results.  An  example 
for  R=1000  m,  Ex=ey=  1  mrad  is  given  in  Figure  5. 


Figure  5  Aimpoint  and  dispersion  generated  shot  lines;  R=1000  m,  e=0.001  rad 

This  kind  of  approach  enables  the  simulation  of  gun  fire  engagements,  e.g.  a  Goalkeeper  gun 
system  defending  a  ship  against  a  sea  skimming  missile. 

The  frequently  used  option  is  a  grid  of  shot  lines  simulating  the  detonation  of  a  warhead,  which 
is  assumed  to  generate  numerous  fragments  emanating  from  the  centre  of  the  warhead. 

Also  the  possible  direct  hit  of  the  intercepting  missile  with  the  target  is  included.  To  this  end  the 
user  must  specify  several  missile  representing  points,  which  are  being  assumed  to  generate 
“direct-hit”  shot  lines,  see  Figure  6. 
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warhead 


Figure  6  Warhead  detonation  and  direct  hit  shot  line  simulation 

This  kind  of  engagement  is  typical  of  an  air  defence  missile,  e.g.  Sea  Sparrow,  defending  a  ship 
against  an  incoming  sea  skimming  missile. 

3.2  Projectile  trajectories 

When  considering  shot  lines  simulating  projectile  trajectories,  the  target  velocity  may  have  a  se¬ 
rious  effect  on  the  impact  condition.  In  general,  it  is  assumed  that  a  projectile  follows  its 
(ballistic)  flight  to  the  target  without  motion  around  its  flight  path,  i.e.  yaw  and  spin  are  ignored. 
However,  at  impact  on  the  target,  flying  with  velocity  Vt,  the  penetration  process  is  changed  sig¬ 
nificantly  due  to  the  velocity  vector  of  the  target  relative  to  the  projectile  velocity  vector,  see 
Figure  7. 


Figure  7  Influence  of  target  velocity 

Given  the  initial  conditions  of  flight,  i.e.  azimuth  (a)  and  elevation  (tir)  angles,  the  original  pro¬ 
jectile  velocity  Vp  is  changed  into  the  relative  velocity  Vpt  by: 


TNO-report 


PML  1994-A33 


Page 

15 


Vpt=Yp+V,  (4) 

The  projectile  now  exhibits  an  angle  of  incidence  to  the  relative  flight  path,  indicated  by  the  dy¬ 
namic  angle  of  incidence  (p,  given  as; 


coscp  = 


(5) 


In  a  first  approximation,  the  shape  r(x)  of  the  projectile  is  assumed  to  be  parabolic  (insert  in 
Figure  8),  i.e.: 

D  rx~  ^ 

r(x)=y  -w  hence:  Amax=  J2r(x)dx  (maximum  presented  area) 


where  r(x)  is  the  local  projectile  radius,  D  the  maximum  diameter  and  L  the  projectile  length. 

The  maximum  presented  area  (A^ax)  is  then  found  as;  A^ax  =  jLD  (viewed  from  the  side), 
whereas  the  minimum  presented  area  is:  Amin=  -  (viewed  from  the  front).  The  ratio  between 
A  and  A^in  is  then  approximated  by: 


A  8  L 

A-—  =  cos(p  -P  3^  Q  sin(p 


(6) 


When  dealing  with  fast  flying  aircraft  or  missiles  and  high-velocity  projectiles,  the  dynamic 
angle  of  incidence  (9)  may  become  as  large  as  20°,  see  Figure  8. 
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Figure  8  Influence  of  dynamic  impact  angle  on  the  presented  area  of  a  projectile 


Hence,  instead  of  a  circular  hole,  the  projectile  will  create  an  oblong  perforation  that  signifi¬ 
cantly  reduces  its  penetration  capability. 
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3.3  Fragment  trajectories 

The  basic  shot  line  data,  such  as  mass,  velocity  and  shape  of  the  fragments  are  usually  obtained 
from  experiments.  Given  the  warhead  velocity  vector  Vm,  fragment  (static)  velocity  vector  Vf 
and  target  velocity  vector  Vj,  the  relative  velocity  vector  of  the  fragment  when  hitting  the  target, 
supposing  that  air  drag  is  not  incorporated,  is  equal  to: 

Yte=Yf  +  y„  +  v,  (7) 

Hence,  the  fragments  will  strike  the  target  with  relatively  high  velocities,  and  may  have  an  arbi¬ 
trary  orientation  at  impact  due  to  tumbling. 

3.4  Influence  of  air  drag 

To  account  for  proper  aerodynamic  drag  during  its  flight  towards  the  target,  each  fragment  is  as¬ 
sumed  to  have  a  specific  average  shape,  expressed  by  its  non-dimensional  average  shape  number 
[5],  [6],  [7]: 


where  Ap  is  the  average  presented  area  and  Vol.  is  the  volume  of  the  fragment. 
Dehn  [6]  gives  an  expression  for  the  average  presented  area  as: 

where  is  the  wetted  surface. 

For  a  number  of  shapes,  the  average  shape  number  can  be  found  as: 

_  Jl 

for  a  sphere :  =  1.21  (independent  of  size) 


(8) 


(9) 


(10) 


for  a  beam : 


S„=l(ab)‘'^  1+-  + 


where  a  =  L/t,  b=w/t,  for  a  beam  with  dimensions  Lxwxt  (L>w>t). 
For  a  cube,  where  all  dimensions  are  equal,  it  yields: 


Sj^  =  1 .5  (independent  of  size) 


For  a  cylinder,  it  can  be  shown  that: 


(11) 

(12) 

(13) 


where  a  =  L/D  of  the  cylinder. 
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When  air  drag  must  be  incorporated,  the  absolute  velocity  of  the  fragment  (i.e.  Ym+Yf)  must  be 
used  to  determine  certain  air  drag  parameters,  e.g.  drag  coefficient,  dynamic  pressure  etc. 
Suppose  the  range  to  the  target  is  R,  and  the  initial  absolute  fragment  velocity  is  Yfo-  The  initial 
relative  velocity  to  the  target  is  then  (Figure  9): 

Yfto=Yf„  +  v^  (14) 

If  the  fragment  is  retarded  by  air  drag,  the  velocity  at  a  range  R  from  the  origin  can  be  estimated 
for  small  velocity  changes,  hence  a  constant  drag  coefficient,  as  [7]: 

Yfi=YfoF(R)  (15) 


or,  when  applied  to  the  relative  fragment-target  velocity: 


Yf,i=V„„.F(R)+V,[l-F(R)] 


where  the  correction  term  F(R)  has  been  derived  in  [7]  as: 


F(R)  =  exp 


-CD-ft 


m 


(16) 


(17) 


where  Pa  is  the  density  of  air,  m  is  the  fragment  mass  and  Cd  is  the  drag  coefficient.  The  effect 
of  burst  altitude  is  included  by  virtue  of  changes  in  air  density.  From  these  velocities,  the  average 
velocity  (i.e.  harmonic  mean)  can  be  derived: 


V  = _ 

— fta  1 


Yfto^Yfti 


—  =  Yfa+Yt 


(18) 


where  Yfa  is  the  average  relative  fragment  ejection  velocity  vector  (Figure  9). 
The  difference  in  time-of-arrival  can  be  approximated  as: 


At  = 


(19) 


Hence,  the  (rearward)  shift  in  impact  due  to  the  retardation  (equal  to  the  difference  in  times  of 
arrival)  is  equal  to: 


As  =  Yt^t  =  RYt 


(20) 


On  average,  the  fragment  originally  ejected  at  an  angle  of: 


-1 


(Yfio-Yi) 

|Yf,„||Y.| 


Xo  =  COS 


(21) 
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is  now  ejected  with  an  average  velocity  Vfja  at  an  average  angle  of; 


Xa  =  cos 


Yf,.||Y,|_ 


(22) 


It  must  be  noted  that  air  drag  becomes  important  when  the  warhead  detonates  long  distances 
from  the  target,  and  hence  retardation  of  fragments  affects  the  flight  trajectories  and  will  cause  a 
shift  of  impact  points  that  may  not  be  ignored,  see  Figure  9. 


Xn^a  Vt 


Figure  9  Ejfects  of  air  drag  on  fragment  trajectories 

In  MISDAC,  it  is  simply  assumed  that  for  the  velocities  involved,  a  constant  Cd  may  be  taken, 
which  is  linearly  interpolated  for  a  given  range  of  Mach  numbers.  The  Mach  number  M  is  de¬ 
fined  as: 

M  =  J^  (23) 

where  c^  is  the  speed  of  sound  of  the  ambient  air.  For  air,  the  speed  of  sound  can  be  expressed  as 

[8]: 


where  7=  1.4  is  the  ratio  of  specific  heats,  Ra=  8314.32  Joule/K.kmol  is  the  absolute  gas  con¬ 
stant,  Mo=  28.9644  kg/kmol  is  the  molecular  weight  of  air  and  T  =  288. 1 5  K  is  the  temperature. 
At  sea  level  it  yields:  Ca  =  340.3  m/s. 

In  [7]  an  equation  is  listed  including  the  average  shape  numbers  for  a  fragment  and  sphere  as: 
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where  Sp 
numbers: 


sphere 


=  1.209  and  Sn 


fragment 


=  1.856^  Table  1  lists  some  Cp  values  for  a  range  of  Mach 


Table  I  Cq  values  for  a  sphere  and  a  (natural)  fragment 


0.0 

0.7 

1.1 

5.0 

9.0 

^^fragment 

1.06 

1.09 

1.40 

1.20 

1.12 

^^sphere 

0.47 

0.57 

1.00 

0.91 

0.91 

The  altitude  of  the  weapon  burst  is  included  in  the  ambient  density  Pa,  which  is  calculated  ac¬ 
cording  to  the  US  standard  atmosphere  [8],  Given  the  velocity  brackets,  a  look  at  Table  1  pro¬ 
vides  the  necessary  data  for  the  drag  coefficient  to  be  used  in  equations  (25),  (16)  and  (17)  to 
calculate  the  impact  velocity  at  the  target.  To  avoid  unnecessary  computations,  a  velocity  deficit 
of  5  pet  is  taken  as  criterion,  or  referring  to  eqn.  (17):  F(R)  =  0.95.  Considering  a  general  frag¬ 
ment  with  moderate  velocity  and  drag  coefficient,  this  results  in  the  criterion  for  range,  beyond 
which  air  drag  calculations  must  be  included  as: 


R 


drag 


>25 


m 


/3 


(m) 


(26) 


Hence,  for  a  fragment  of  m=0.01  kg  and  Sn=l  .9  this  range  becomes  approximately  3  m,  but  for  a 
large  fragment  (m=0.08  kg,  Sn=1.9)  this  range  becomes  approximately  5.6  m. 

Effects  of  gravity  are  excluded,  since  the  trajectories  considered  are  relatively  short  for  gravity  to 
become  important. 


3.5  Influence  of  fragment  tumbling 

To  bring  the  tumbling  of  a  fragment  into  account,  the  presented  area  of  the  fragment  at  the  mo¬ 
ment  of  impact  (Ap)  is  assumed  to  vary  slightly  from  its  predetermined  value  between  A^jn  and 
A-max- 

A  procedure  to  estimate  the  presented  area  for  a  beam-shaped  fragment  (at  impact),  viewed  from 
a  surrounding  sphere  with  equally  distributed  viewpoints,  resulting  in  a  non-homogeneous  area 
distribution  is  proposed  in  [9],  see  Figure  10. 


9  1 

In  ref.  [7]  the  US-based  shape/ac/or  C  [ft  grain  /lb]  is  used,  the  conversion  being:  8^^=  3.25C. 


1 
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A/A 

average 


Figure  10  Presented  area  ratio  histogram  (  Sfi=l  .856  ,  based  on  /  000  samples) 


The  presented  area  at  impact  defines  the  impact  shape  number: 


S 


n 


A 

Vol 


P _ 

2/3 


(27) 


An  example  for  three  different  shape  numbers,  i.e.  for  a  sphere,  a  beam  and  a  cube  (with  equal 
masses),  is  given  in  Figure  1 1 . 
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Figure  1 1  Ejfects  of  tumbling  on  impact  shape  number  distribution 


This  approach  demonstrates  the  possibility  that  from  a  class  of  fragments  having  an  average 
shape  number,  which  makes  it  impossible  to  penetrate  a  certain  target  plate,  one  or  more 
impacting  fragments  will  still  be  able  to  penetrate  this  plate  due  to  the  (randomly)  reduced  im¬ 
pact  shape  number. 
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In  MISDAC,  each  ‘‘shot  line”  representing  a  fragment  trajectory  has  a  random  shape  number  as¬ 
signed  to  it  corresponding  with  the  random  actual  presented  area  at  (first)  impact.  Each  fragment 
is  assumed  to  have  a  beam-like  shape,  which  is  expressed  by  the  average  shape  number: 


S„=^(ab)''^  1  +  --P- 


where  a=L/t  and  b=w/t. 

A  procedure  to  estimate  the  values  of  a  and  b  for  an  arbitrary  (natural)  fragment  with  a  given 
shape  number  Sp  is  [9]: 


a  =  l  +  13.183(Sn  -1.5) 

b  =  l+  4.700(Sn -1.5) 


(S„>1.5) 


The  “optimal”  values  for  a  and  b  to  match  the  beam  with  the  given  shape  number  can  be  found, 
using  the  “scaling”  expressions,  assuming  a  constant  ratio  of  a/b: 

a  =  aoX  and  b  =  boX  (30) 

where  ao  and  bp  are  the  original  constant  values  of  a  and  b.  A  Newton-Raphson  iteration  is  now 
used  to  find  x,  which  estimates  the  roots  of  an  equation  according  to  the  iterative  relation: 
f(x) 

Xn+l=Xn‘f^  ri=f . m  iterations  (31) 

where: 


f(x)  =  Sn-{(aobox2)'^3  ^ 

aoX  boX 

r  n  (32) 

dx  [x  a„x2  baX^j 

This  allows  a  quick  assessment  of  the  appropriate  values  for  the  equivalent  beam-shaped  frag¬ 
ment.  It  should  be  noted,  that  inserting  the  value  Sn  =1.5,  which  is  the  value  for  a  cube  where 
a/b=l,  gives  a  trivial  solution.  Hence,  the  solutions  found  are  valid  for  Sn  >  1.5. 

For  further  penetration  in  the  target,  the  value  of  the  shape  number  is  assumed  constant  for  this 
trajectory  since  the  effect  of  rotation  velocity  is  neglected.  If  it  assumed  that  the  fragment  rotates 
(tumbles)  at  w=3000  rpm^  then  for  a  lateral  velocity  of  V=1000  m/s,  the  distance  for  a  complete 
rotation  would  be  As=V/(o=20  m,  which  is  much  greater  than  the  usual  distances  encountered  in 
(aircraft)  target  (component)  dimensions,  apart  from  the  maximum  distance  a  fragment  can 
penetrate  into  this  target. 


This  value  is  considered  extremely  fast  by  ballistic  experts  at  TNO-PML. 
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3.6  Impact  and  penetration 

When  the  target  is  hit  by  the  fragment  (or  projectile),  the  penetration  capacity  must  be  assessed. 
In  general,  for  steel  fragments  hitting  a  target  plate,  the  general  THOR  equation  [10]  is  used,  ex¬ 
pressed  by: 

Vr  =  V-10‘^l[tApJ'^m‘'3sec6‘'4vC5  (33^ 

where  is  the  residual  velocity,  t  is  the  target  plate  thickness^,  Ap  is  the  presented  area  at  im¬ 
pact,  0  is  the  normal  impact  angle  (sec0=l/cos9)  and  m  and  V  are  the  mass  and  velocity  of  the 
fragment  at  impact.  The  THOR  constants  c  j  -  c ^  are  based  upon  experiments.  Plotted  on  a  log- 
log  scale  eqn.  (33)  provides  linear  terms,  which  can  easily  be  used  for  fitting  the  constants  with 
experimental  data  (e.g.  multiple  regression).  The  variables  used  and  the  shape  of  eqn.  (33)  does 
not  have  a  physical  background,  except  for  their  importance  in  ballistic  experiments.  It  must  be 
noted  that  the  use  of  the  constants  are  valid  only  within  the  limits  of  the  experiments,  and  ex¬ 
trapolation  to  higher  velocities,  for  example,  will  introduce  errors.  In  fact,  when  the  THOR 
equation  is  used  for  a  higher  velocity  regime  than  the  constants  have  been  derived  for,  new  ex¬ 
periments  are  required  or  erroneous  results  must  be  accepted. 

The  possibility  of  ricochet  is  included  by  limiting  the  impact  angle  -7O°<0<7O°  [11].  Shot  lines 
impacting  outside  these  boundaries  are  rejected  (Figure  12). 


Figure  12 


Using  the  impact  shape  number  (eqn.  27),  eqn.  (33)  can  be  rewritten  as: 

Vr  =  V-10‘'l[tSn  sec0‘^"^V‘^5  (34) 

where  Pni=7850  kg/m  is  the  density  of  steel,  since  THOR  is  valid  for  steel  fragments.  The 
maximum  penetrable  plate  thickness  -  t^ax  ■  is  found  by  setting  Vr=0,  which  yields: 


“^max 


^  jQ-ci  /c2 


.2/3 


m 


-2/3-C3/C2  cec0~‘^4/c2' 


(35) 


This  thickness  is  equal  to  the  maximum  distance  (d^ax)  the  fragment  can  penetrate  into  this  tar¬ 
get  plate,  while  in  general:  d  =  tsec0,  where  0  is  the  normal  impact  angle  and  d  is  the  line-of- 
sight  distance.  From  eqn.  (35)  it  can  be  seen  that  tmax  is  inversely  proportional  with  the  shape 
number  Sn-  To  illustrate  the  importance  of  the  shape  number.  Figure  13  shows  the  effect  of  the 


2 


Not  to  be  confused  with  penetration  distance,  or  line-of-sight  distance. 
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(impact)  shape  number  on  the  maximum  penetration  capability  for  a  number  of  impact  velocities 
of  a  fragment  mass  of  25  gram  impacting  on  an  aluminium  plate  at  normal  obliquity. 


3.7  Blast  effects 

The  effects  of  blast  are  hard  to  implement  in  a  shot  line  oriented  vulnerability  code.  A  major 
variable  in  blast  related  analyses  is  the  miss  distance  between  the  warhead  burst  point  and  the 
target.  In  MISDAC  the  blast  vulnerable  component  is  modelled  as  a  line  element,  for  which  the 
miss  distance  can  be  easily  derived.  In  parameter  form  the  equation  describing  the  position  and 
orientation  of  the  line  rn  yields: 

x  =  M  +  >.m  (36) 

where  ^  is  a  parameter,  M  is  a  fixed  point  on  the  line  element  (usually  the  centre  position)  and  B 
is  the  warhead  burst  point  (Figure  14). 
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Figure  14  Definition  of  the  miss  distance 

The  line  BX=  X-B  is  normal  to  m  if  the  scalar  vector  product  equals  zero,  or: 

(^•m)  =  0  (37) 


or: 


r  ^ 

< 

M  -1-  A.rn  -  B 

•  m 

X  '  J 

(38) 


from  which  A,  can  be  solved  as: 

{(B-M)-m} 

(m»m) 

The  length  of  the  vector  BX*  is  then  equal  to  the  miss  distance  MD  to  be  found  as: 
MD  =  |M  +  A.*m-B| 


(39) 


(40) 


The  miss  distance  is  a  major  variable  which  defines  the  encounter  geometry  (see  also  Figure  2). 
When  a  warhead  detonation  occurs,  the  blast  wave  generated  is  a  spherical  expanding  wave, 
which  can  be  characterised  by  its  peak  (side-on)  overpressure  ps  and  positive  duration  tp,  see 
Figure  15.  In  general,  the  relation  between  peak  overpressure  and  time  is  expressed  by  the  modi¬ 
fied  Friedlander  equation  as: 


P(t)  =  Ps 


exp 

^P- 

(41) 


where  B  is  a  slowly  varying  time  factor,  which  is  considered  constant  for  this  situation.  A  similar 
equation  is  used  to  describe  the  relation  of  the  flow  velocity  or  particle  velocity  as  a  function  of 
time.  For  t  >  tp,  the  pressure-time  curve  exhibits  a  negative  phase,  or  suction  phase,  where  the 
motion  of  (air)  particles  is  reversed  and  is  directed  to  the  explosion  source. 
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Figure  15  Blast  wave  characteristics 


The  motion  of  particles,  denoted  by  their  (flow)  velocity  u,  will  also  create  a  dynamic  pressure 
when  their  flow  is  distorted.  If  the  flow  would  be  isentropically  brought  to  rest  a  dynamic  pres¬ 
sure  q=const.pu2  is  created,  where  p  is  the  flow  density.  In  a  first  approximation,  the  dynamic 
pressure  originates  from  the  one-dimensional  Euler  equation  of  motion  [12]: 


0U 


du  1  3p 

+  u  —  + - - 

dx  p  d\ 


(42) 


If  the  pressure  wave  would  propagate  without  distortion  and  any  decay,  it  can  be  written: 

3u  _  du 
3t  dx 


(43) 


where  c  is  the  wave  speed  (constant). 

Inserting  this  term  yields,  assuming  a  constant  density  (incompressible): 


which  can  be  integrated  to  yield  the  equation: 


p  -I-  jpu^  +  puc  =  constant 

V - V 

Stagnation 

pressure 

(Bernoulli) 


(45) 


This  equation  differs  substantially  from  the  Bernoulli  equation  [12],  which  omits  the  term  puc, 
but  which  is  valid  for  a  steady  incompressible  flow,  i.e.  3u/3t  =  0.  Basically,  when  the^ highly  un¬ 
steady  blast  flow  field  is  disturbed,  the  (unsteady)  dynamic  pressure  -  indicated  by  -pu^+puc  - 
will  be  much  larger  than  originally  conceived  from  the  Bernoulli  equation.  However,  blast  waves 
generally  are  not  simple  waves  having  a  constant  amplitude  and  speed.  Table  2  lists  a  number  of 
values  for  a  spherical  charge  of  TNT  of  50  kg,  which  shows  how  fast  a  blast  wave  decays. 
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target  at 
intercept 


target  at  morrent 
of  warhead  burst 


Figure  16  Effect  of  target  velocity 
The  shock  wave  velocity  can  be  approximated  as: 
R 


U,=a 


W 


1/3 


(47) 


of  which  the  constants  are  found  using  regression  techniques.  The  average  shock  wave  velocity 
can  then  be  derived  as: 


Us=— ^f^U3(x)dx  (48) 

L-  •'Zc 

where  Z=R/W''^^  is  the  scaled  range. 

Using  the  expression  of  Us,  eqn.  (47),  the  approximation  yields  for  0.053^<1 .59: 

-  2962  (m/s)  (49) 

Z -0.052  (  \  ^  ’ 


The  minimum  range  which  should  be  observed  is  equal  to  the  charge  radius: 
Rc=  0.053W’^^  (m),  at  which  a  “direct  hit”  or  contact  explosion  occurs. 

Based  upon  geometrical  considerations  (Figure  16)  the  distance  at  which  the  blast  strikes  the 
target  can  be  found  as  [14]: 


2 

1 

VjMD 

1 

MD 

f 

U.R; 

-ui| 

.  Ri  . 

(50) 


where  Rj  is  the  initial  distance  at  the  moment  of  intercept.  The  “+”  sign  refers  to  an  approaching 
target  whereas  the  sign  refers  to  a  target  moving  away.  The  average  velocity  U  is  based  on 
the  original  range  Rj.  With  respect  to  Figure  16,  it  can  be  shown  that  the  sign  is  equivalent  to  the 
scalar  vector  product: 

sign  =  -  (BT»Vt)  (>  0  for  an  approaching  target). 


(51) 
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As  an  example,  consider  a  50  kg  charge  of  TNT  detonating  at  an  intercept  range  of_5  m  from  an 
approaching  target  flying  at  300  m/s  at  a  miss  distance  of  3  m.  Then  it  follows  that:  U  =  1935  m/s 
and  Rhit~  4.5  m,  or  a  10%  reduction  in  range.  This  will  give  a  distinct  increase  in  peak 
overpressure,  since  the  scaled  range  R/W*'^^,  is  reduced  from  1.36  to  1.22  m/kg*'^^,  respectively. 
The  original  peak  static  overpressure  will  increase  from  430  to  600  kPa,  which  is  an  increment  of 
almost  40%,  because  of  the  exponential  behaviour  of  the  overpressure -distance  curve. 

Apart  from  the  virtual  “closer”  range,  the  blast  wave  will  hit  a  moving  target,  which  will  have  its 
impact  on  the  dynamic  reflection  process.  In  general,  shock  wave  reflection  relations  can  be  es¬ 
tablished  for  a  stationary  situation,  but  this  changes  dramatically  when  the  target  is  moving.  In 
addition,  the  blast  loading  duration  on  the  target,  defined  by  the  time  the  shock  wave  needs  to 
pass  the  target^  will  be  reduced  conversely  (similar  to  the  Doppler  effect).  Figure  17  illustrates 
these  effects. 

It  is  assumed  that  an  ideal  reflection  takes  place  of  a  shock  wave  against  an  opposed  moving  wall 
with  velocity  -V,  where  Ma  and  Mb  are  the  Mach  numbers  of  the  incident  and  reflected  shock 
waves  in  regions  1  and  2,  respectively  (Figure  17). 

stationary  moving 


Figure  1 7  Reflection  against  a  stationary  (a)  and  moving  wall  (b),(c) 

A  -  incident  shock  wave,  B  -  reflected  shock  wave 

A  simple  geometric  consideration  with  respect  to  the  boundary  condition  at  the  wall  yields  the 
reflection  condition  [12]: 

U3-U2  =  -(u2-ui)-i-V  (52) 


where  V  is  the  wall  velocity  (positive  to  the  right). 

The  Rankine-Hugoniot  equations  are  used  to  express  U2,  C2  and  U3  in  Ma  and  Mb: 


U2-»l 

Cl 

U3  -U2 
C2 


2 

7  +  1 

> 

1 

> 

right  travelling  wave 

> 

-2 

1 

left  travelling  wave 

7  +  1 

Mb-— - 

Mb_ 

(53) 


1 


Called  the  diffraction  loading. 
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Using  the  reflection  condition,  eqn.  (52),  Mb  can  then  be  solved  from  the  resulting  quadrature: 


-Mb^ 


Sl 

C2 


2C2 

dynamic 

reflection 


-1  =  0 


(54) 


from  which  the  reflected  pressure  p3  is  calculated  according  to: 


(55) 


The  reflected  pressure  p3  can  be  put  together  in  the  reflection  coefficient  Rc,  defined  as; 

P2-P1 


(56) 


Using  the  expressions  for  Ma  and  Mb,  which  yield  solutions  for  p2  and  P3,  the  reflection  coeffi¬ 
cient  can  be  found  as  a  function  of  the  initial  shock  strength  (or  Ma).  Omitting  the  details. 
Figure  18  shows  the  reflection  coefficient  for  V=±300  m/s,  which  shows  a  dramatic  difference 
with  “normal”  (i.e.  V=0)  reflection.  When  V  is  positive,  an  expansion  may  occur  immediately 
behind  the  reflected  shock  wave,  which  is  the  case  when  the  flow  velocity  U2  is  smaller  than  the 
wall  speed.  This  phenomenon  has  been  displayed  in  Figure  17  (c).  In  that  case,  the  reflection 
overpressure  decreases  and  even  underpressure  may  develop. 


Figure  18  Dynamic  reflection  coefficient 

It  should  be  noted  that  Rc  increases  at  low  overpressure  ratios  (for  an  opposing  wall),  which  is 
caused  by  the  relatively  large  influence  of  V  compared  with  the  particle  velocity.  Conversely, 
this  coefficient  decreases  for  an  expanding  wall  movement  at  low  overpressure  ratios.  Obviously, 
when  the  wall  speed  is  equal  to  the  shock  wave  velocity  itself,  reflection  does  not  occur  at  all. 
Another  interesting  feature  is  the  diffraction  impulse,  which  is  controlled  by  the  time  the  blast 
waves  needs  to  “pass”  the  target.  In  general,  this  time  can  be  expressed  in  terms  of: 
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“  Us+V  (5"^) 

where  L  is  an  arbitrary  dimension.  The  product  of  reflected  pressure  and  diffraction  period  is 
then: 

Id  =  Retd  (58) 

which  can  be  considered  to  represent  the  diffraction  impulse.  Figure  19  shows  as  a  function  of 
the  incident  shock  strength,  and  clearly  shows  that  the  net  result  will  be  a  dramatic  increase  of 
impulse  compared  with  the  stationary  target  (V=0). 


Figure  19  Reflection  (diffraction)  impulse 

Despite  the  shorter  duration  (similar  to  the  Doppler  effect),  the  reflected  pressure  is  much 
stronger  leading  to  a  significant  increase  in  loading  on  the  (opposed  moving)  target.  Obviously, 
this  effect  will  be  reversed  when  the  target  moves  away  from  the  incoming  shock  wave,  but  this 
effect  requires  complex  computations  including  interaction  of  rarefaction  waves  with  shock 
waves,  and  is  therefore  omitted  in  this  analysis. 

MISDAC  incorporates  a  blast  kill  routine  (see  paragraph  4.4)  by  assuming  that  the  user  provides 
blast-critical  elements  for  the  target.  Each  element  is  a  line  element,  i.e.  a  cylindrical  solid  (with 
vulnerability  code  2)  with  zero  radius,  so  that  only  its  orientation  and  length  provide  the  essential 
information  with  respect  to  warhead  burst  point.  For  example,  the  vulnerability  of  the  fuselage  of 
a  missile  to  blast  is  simply  represented  by  a  line  in  its  centre,  with  a  certain  length.  The  code  then 
determines  the  miss  distance  and  the  range  of  this  element  to  the  detonation  point  (which  is  as¬ 
sumed  to  occur  in  the  centre  of  the  warhead)  and  accounts  for  relative  speeds  to  determine  the 
range  when  hit  by  the  blast.  A  structural  analysis,  however,  is  not  performed,  since  this  requires 
knowledge  of  blast  loading  histories  (a  problem  on  its  own),  structural  layout  (for  example,  a  fi¬ 
nite  element  analysis)  and  dynamic  structural  response  analysis,  which  can  usually  be  performed 
by  advanced  finite  element  codes.  This  combination  is  too  cumbersome  to  apply  in  MISDAC, 
where  use  is  made  of  simple  analyses  requiring  characteristic  properties.  In  general,  when  deal¬ 
ing  with  blast  damage  aspects,  it  is  sufficient  to  employ  global  methods,  such  as  the  damage  ra¬ 
dius,  which  can  be  expressed  by  the  Hopkinson  scaled  range  R/W’^3^  irrespective  of  target  orien¬ 
tation,  size  etc.  Besides,  it  should  be  noted  that  blast  damage  is  usually  a  secondary  damage 
mechanism,  since  fragment  damage  dominates  the  overall  damage  capability  of  fragmenting 


TNO-report 


PML  1994-A33 


Page 

31 


warheads.  Paragraph  4.4  deals  with  the  assumption  put  forward  in  the  blast  damage  analysis  such 
as  used  in  MISDAC. 


4  DAMAGE  AND  KILL  CRITERIA 


A  major  item  in  missile  vulnerability  assessment  is  the  assignment  of  the  appropriate  damage 
criteria  to  the  respective  components  considered.  In  MISDAC,  four  options  exist  to  define 
specified  damage  criteria,  viz.: 

•  Penetration  damage  (fragmenfprojectile) 

In  general,  this  mechanism  pertains  to  components  having  specific  functions  such  as 
power  transfer  or  to  components  that  have  not  been  described  in  detail,  but  in  which  cer¬ 
tain,  “processes”  take  place,  such  as  electronics,  gear  box  etc.  For  this  kind  of  damage  cri¬ 
terion,  two  options  are: 

•  complete  penetration  (perforation),  e.g.  holes  in  plates  (fuel  tanks)  etc; 

•  minimum  penetration  distance,  e.g.  plate  thickness  to  be  penetrated.  For  this 
mechanism,  the  tso  criterion  is  used  (see  later). 

•  Impact/Shock  damage 

This  mechanism  occurs  when  a  slug  hits  the  component.  The  impact  (kinetic)  energy  may 
cause  shock  damage,  impact  cratering  and  initiate  other  processes,  such  as  shock  initiation 
of  explosives.  For  this  mechanism,  the  critical  kinetic  energy  criterion  (E50)  is  used. 

•  K(ill)-factor  curves  (direct  component  kill) 


•  Blast  (shock  wave)  damage 

This  mechanism  applies  when  the  warhead  contains  a  large  explosive  charge  (e.g.  nuclear 
weapons,  Fuel-Air  explosives,  large  conventional  weapons)  whose  primary  effect  (or  in 
addition  to  fragment  effects)  is  the  blast  wave  striking  the  missile  in  flight.  In  general,  the 
kill  by  blast  is  the  structural  kill.  For  this  mechanism,  the  critical  radius  to  the  charge  (R50) 
is  used  as  criterion. 

Several  options  stem  from  the  code  TARVAC,  but  have  been  adapted  for  the  purpose  of  missile 
damage  assessments. 

Basically,  the  probability  of  kill  given  a  hit  on  a  component  (P^h)  by  the  threat  is  defined  by: 

Pkh=Pdh-Pkd  (59) 


where: 

Pdh  probability  of  damage  given  a  hit 
Pkd  probability  of  kill  given  damage 

The  first  variable  depicts  the  physical  damage  to  a  component,  whereas  the  second  variable  more 
or  less  depicts  the  functional  degradation  of  a  component.  Its  combination,  P^h^  does  not  reflect 
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the  probability  of  kill  of  the  complete  missile,  however,  since  this  incorporates  a  systems  kill  in 
conjunction  with  residual  flight  capability  to  hit  the  target.  Hence,  eqn.  (59)  reflects  the  proba¬ 
bility  of  kill  for  the  internal  components  of  a  missile  target  (or  any  target),  which  will  have  fur¬ 
ther  consequences  for  the  missile's  residual  flight  trajectory. 

The  probability  of  damage  (of  a  component)  is  estimated  given  the  damage  mechanism  (e.g. 
penetration)  by  a  relation  including  the  estimated  critical  value  for  that  particular  criterion,  which 
is  approximated  by: 


X/X50  +exp[-|(x/x5o-l)] 


where  X50  is  the  estimated  critical  value  indicating  a  50%  probability  when  the  criterion  (i.e. 
x=X5o)  is  met.  This  relationship  is  quite  close  to  the  normal  cumulative  probability  distribution. 
The  probability  function  accounts  for  the  uncertainties  in  the  assessment  since  the  exact  location 
of  the  hit  and  the  direction  the  penetrator  is  going  are  averaged  over  the  whole  component,  which 
virtually  may  have  more  than  one  vulnerable  area.  The  variable  t,  defines  the  deviation  around 
X50,  or  in  fact  the  “accuracy”  (or  error)  in  the  estimated  critical  value.  Figure  20  shows  the  proba¬ 
bility  curve  for  a  number  of  values  of 


Figure  20  Damage  probability  curves  for  some  ^  values 

It  can  be  shown  that  |  can  be  related  to  the  standard  deviation  (a)  of  a  normal  distribution  as 
[15]: 

^  )^-ln(a  +  l) 


a 


(61) 
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4.1  Penetration  damage  criteria 

For  a  penetrator,  the  penetration  criterion  is  largely  based  on  experimental  data  listed  by  the  Thor 
penetration  relation,  valid  for  steel  fragments,  where  the  maximum  penetration  distance,  see  tmax 
from  eqn.  (35),  is  compared  with  the  estimated  critical  value  tsg. 

The  probability  of  kill  given  damage,  Pkd>  can  be  estimated  utilising  the  Rayleigh  (or  Weibull 
(2,6))  distribution  [16]: 

Pkd  =l-exp[-(y/6f  ]  (62) 

where  6  =  yso/V  ln(2)  is  the  scale  factor  (with  median  yso),  or: 

Pkd  =  l-exp -ln(2)(y/y5of  (63) 

which  gives  a  50%  probability  when  y=y5o-  This  curve  very  much  resembles  the  damage  proba¬ 
bility  function  when  ^=2  (equivalent  with  a  normal  distribution  with  0=0.6).  The  probability 
function  accounts  for  the  uncertainties  in  the  assessment  of  the  (functional)  kill  of  the  compo¬ 
nent,  which  may  depend  on  the  location  of  the  damage  and/or  the  extent  of  damage  inflicted.  In 
eqn.  (63)  y/yso  is,  for  example,  the  effect  of  the  penetration  hole  size  (A/A50)  on  the  actual  kill  of 
the  component  given  damage  by  a  penetrator.  For  example,  when  A=0  while  Pdh=l  then  still 
Pkd=0,  because  a  penetrator  with  an  infinitely  small  penetration  hole  area  cannot  produce  a  kill. 
It  may  be  noted  that  eqn.  (60)  may  be  used  here  also,  but  then  the  variable  ^  must  be  specified, 
which  introduces  an  extra  variable. 

4.2  Impact/Shock  damage  criteria 

This  mechanism  is  quite  similar  to  the  penetration  process,  where  a  fragment/projectile  or  a  par¬ 
ticle  hits  a  component  with  a  certain  mass  and  velocity.  Whether  or  not  penetration  or  perforation 
is  achieved,  the  impact  kinetic  energy  may  cause  internal  shock  damage  (mechanical  accelera¬ 
tions),  cratering  or  some  other  deformation.  As  the  x^q  criterion,  the  impact  kinetic  energy  = 
l/2mV^  is  used,  in  addition  to  the  probability  equation  (60).  Similarly  to  the  previous  mecha¬ 
nism,  the  probability  of  kill  given  damage  is  assessed  with  the  critical  area  A^q  as  the  (y^Q)  cri¬ 
terion,  since  it  is  assumed  that  energy  density  (in  fact  E,-„/A^^)  is  important  in  controlling  the 
kill,  see  [20]. 

4.3  Direct  kill  criteria 

A  direct  attempt  to  assess  the  kill  of  a  component  given  a  hit  is  the  employment  of  the  pre-de- 
termined  vulnerability  curves,  indicated  by  a  so-called  K(ill)-factor  [3],  The  larger  the  value  of 
K,  the  more  vulnerable  the  component.  At  present,  similar  to  TARVAC,  MISDAC  allows  for  up 
to  nine  levels,  expressed  by  the  variable  log(mV  )  (Figure  21).  The  curves  are  categorised  from 
“practically  invulnerable”  (K=l)  to  “highly  vulnerable”  (K=9)  given  a  (single)  hit.  Selection  of  a 
specific  value  for  a  component’s  vulnerability  should  be  preferably  based  on  experimental  data. 
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log{mv®) 


Figure  21  Pre-determined  kill  probability  curves  (K-factor) 

In  general,  the  kill  probability  given  a  hit  using  the  K-factor  is  assumed  by: 

Pkhc  ~ 

1  -I-  exp 

and  specific  values  for  the  constants  a,  b  and  c  apply  for  various  K-values,  see  Table  3. 

This  expression  allows  for  the  possibility  that  a  component  cannot  be  killed  by  a  (single)  frag¬ 
ment,  e.g.  such  as  may  happen  with  a  large  electronic  unit  being  hit  by  a  relatively  small  frag¬ 
ment. 

Table  3  Values  of  vulnerability  constants 


K 

a 

b 

c 

1 

0.3 

9.0 

9.0 

2 

0.4 

8.5 

8.0 

3 

0.5 

8.5 

7.0 

4 

0.6 

9.0 

6.0 

5 

0.7 

8.0 

5.5 

6 

0.8 

6.0 

4.9 

7 

0.9 

6.0 

4.3 

8 

1.0 

6.0 

3.5 

9 

1.0  . 

7.0 

2.7 

Some  experimental  data  exist  that  nicely  correspond  to  a  specific  K-value,  but  at  the  moment  the 
basic  principle(s)  could  not  be  verified,  other  than  being  the  result(s)  of  extensive  data  analyses. 
One  analysis  indicates  that  the  variable  mV^  -  being  rewritten  as  VmO.33  .  corresponds  quite  well 
with  associated  Thor  constants  for  V  and  m  (1.0  and  0.333,  respectively)  for  some  materials  used 
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in  helicopter  vulnerability  studies,  in  which  penetration  criteria  (tmax)  based  on  Thor  penetration 
relations  were  used. 

In  general,  these  curves  may  be  applied  for  the  kill  assessment  of  complicated  (and  possibly 
large)  components,  where  a  single  fragment  causes  damage,  but  not  sufficient  to  produce  a  kill. 
Also,  quick  analyses  can  be  performed  when  a  target  is  composed  of  a  large  number  of  compo¬ 
nents  (e.g.  a  ship). 

The  actual  kill  of  a  component  given  multiple  hits  by  fragments  can  be  assessed  by  the  combina¬ 
tion  of  their  individual  kill  probabilities: 

N 

Pkc=i-n(i-Pkj)  (65) 

j=i 


4.4  Blast  damage  criteria 


As  explained  earlier,  a  structural  analysis  cannot  be  performed  in  a  shotline  oriented  code  such  as 
MISDAC.  However,  it  is  sufficient  to  use  global  methods  to  account  for  structural  response  by 
employing  the  so-called  blast  damage  radius,  or  critical  blast  radius.  Basically,  the  target  re¬ 
sponse  is  assumed  to  follow  the  P-I  concept,  where  P  stands  for  (maximum)  applied  load  and  I 
stands  for  (maximum)  applied  impulse.  Any  mechanical  structure  responds  in  a  special  manner 
to  an  applied  load,  but  stays  within  two  limiting  realms,  namely  the  quasi-static  loading  realm 
where  the  structural  response  depends  only  on  P  and  the  impulse  loading  realm  where  the  struc¬ 
tural  response  depends  only  on  the  impulse  I.  For  any  response  between  these  boundaries,  called 
the  dynamic  realm,  the  structural  response  depends  both  on  P  and  I.  It  can  then  be  shown 
[14], [17]  that  blast  damage  can  be  related  to  the  main  blast  damaging  parameters,  viz.  the  range 
to  the  warhead  (R)  and  the  warhead  size  (W). 

For  (external)  blast,  the  critical  radius  Rent  can  be  expressed  by  the  Westine  (R-W)  vulnerability 
relation  [17]: 


^crit 


AW 


1/3 


nl/6 


1  H - 1 - TT 

w 


(66) 


where  A,  B  and  C  are  component  dependent  constants^  and  W  is  the  equivalent  (spherical)  free- 
air  TNT  charge  weight.  With  the  three  constants  given,  it  represents  the  entire  range  of  loads  of 
the  blast  wave  on  the  structure.  The  equation  covers  the  field  from  quasi-static  (or  pure  pressure) 
loading  (B=C=0)  to  combined  pressure-impulse  loading  (C=0)  up  to  purely  impulsive  loading 
(B=0).  It  must  be  noted  that  this  equation  may  be  regarded  as  a  general  vulnerability  equation, 
since  it  implies  target  size,  orientation  and  damage  or  failure  mode,  i.e.  the  extent  and  type  of 
damage.  In  general,  selecting  specific  values  for  the  constants  A,  B  and  C  implies  one  particular 
vulnerability  curve  in  the  R-W  plane,  but  the  target  itself  may  be  composed  of  several  curves. 


Similar  to  the  Thor  constants  used  in  penetration  analyses,  the  author  suggests  these  constants  to  be 
called  the  Westine  constants. 


TNO-report 


PML1994-A33 


Page 

36 


where  the  envelope  determines  its  ultimate  vulnerability.  Also,  no  distinction  is  made  between 
blast  generated  by  internal  or  external  explosions,  or  the  effects  of  protection  of  a  target  by 
neighbouring  components  (shielding),  which  in  fact  obstructs  the  line-of-sight  from  the  burst 
point  to  the  target.  However,  incorporating  these  effects  would  require  a  dynamic  target  model, 
which  is  presently  not  at  stake. 

When  dealing  with  high-velocity  weapon-target  interactions,  such  as  in  the  case  of  a  SAM  war¬ 
head  and  sea  skimming  missiles,  the  (R,W)  relation  can  be  compensated  for  velocity  effects  [14], 
[18],  [19].  Since  a  blast-critical  component  is  modelled  as  a  simple  line  element  (i.e.  cylinder 
with  r=0),  the  critical  radius  R^rit  is  identical  with  the  critical  miss  distance  M50.  The  probability 
of  damage  given  a  hit  (Pdh)  is  estimated  using  an  approximation  similar  to  eqn.  (60): 

^dh  =  1 - F— - (67) 

x-l-exp[-^(x-l)] 

where  x  =  MD/M50  and  where  ^  may  be  selected  as  2.0  (by  default).  Figure  22  shows  the  dam¬ 
age  probability  function  for  a  number  of  ^  values. 


Figure  22  (Blast)  damage  probability  functions  for  some  ^  values 

The  probability  of  kill  given  damage  can  be  derived  taking  the  position  of  the  burst  point  relative 
to  the  centre  of  the  line  element  into  account,  which  in  fact  accounts  for  the  presented  area.  The 
probability  of  kill  given  damage,  Pk(i>  is  estimated  utilising  the  Rayleigh  probability  density  dis¬ 
tribution: 


Pkd  -  exp 


i(R/R 


■ln(2)(R/R5of]=[±f  (68) 

where  R50  =  V  Mio+(5L)2  ,  L  is  the  length  of  the  line  element  and  R  is  the  distance  from  the 
centre  of  the  line  element  to  the  burst  point.  This  approach  assumes  that  the  kill  given  damage  by 
blast  is  related  -  in  some  way  -  to  the  way  it  interacts  with  the  target,  e.g.  face-on  or  side-on  with 
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respect  to  the  general  target  dimensions.  For  example,  a  long  rod  or  cylinder  is  quite  sensitive  to 
the  exact  location  of  the  charge,  while  a  sphere  (or  box)  will  be  less  sensitive.  When  dealing  with 
very  large  planes,  e.g.  wing  panels,  it  may  be  appropriate  to  construct  this  component  with 
several  line  elements.  However,  it  should  be  noted  that  characteristic  features  should  be  used. 
Hence,  when  dealing  with  wings,  for  example,  the  main  wing  spar(s)  should  be  modelled  since 
this  structure  carries  the  main  load. 
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5  Systems  damage  evaluation 


The  final  assessment  of  probability  of  kill  of  the  target  can  be  performed  when  it  is  known 
whether  the  (target)  missile  will  hit  its  intended  (ship)  target.  To  this  end,  an  estimate  must  be 
given  of  the  damage  inflicted  to  mission-sensitive  systems,  such  as  autopilot,  aerodynamic  con¬ 
trol  systems  etc.  Obviously,  this  is  the  most  difficult  part  of  the  program,  since  an  evaluation 
must  be  performed  of  the  functional  degradation  of  relevant  systems,  and  its  impact  on 
continuing  flight  performance.  This  requires  thorough  knowledge  of  interior  systems,  and  their 
functional  relations  with  other  systems  on  board.  For  example,  when  simulations  indicate  that  the 
guidance  unit  will  be  hit  several  times  by  a  number  of  fragments  from  a  specific  warhead,  one 
should  establish  the  effects  it  could  impose  on  the  flight  computer  output,  i.e.  steering  signals  to 
the  rudder  actuators.  What  will  happen  when  signals  are  lost,  e.g.  caused  by  severance  of  the  sig¬ 
nalling  cables  to  the  rudder  actuators?  Also,  in  some  instances  some  systems  may  become  more 
important,  and  hence  are  more  susceptible  to  failure.  If,  for  instance,  the  warhead  would  be 
hammered  such  that  a  low  order  detonation  (or  even  high  order)  is  likely,  then  large  pieces  or 
debris  of  the  missile  will  follow  a  ballistic  trajectory  and  may  inflict  damage  to  the  ship.  To  illus¬ 
trate  the  difficulties,  suppose  that  the  sea  skimming  missile  performs  a  dog-leg  manoeuvre  and  is 
hit  during  the  first  leg,  just  before  turning  to  the  ship.  If  radar  contact  is  present,  but  when  con¬ 
trols  fail  and  (supposedly)  take  a  fixed  position,  then  the  missile  would  miss  the  ship  entirely.  If 
the  missile  is  hit  in  the  second  leg  it  is  already  on  course  and  may  possibly  hit  the  ship. 

This  part  will  be  addressed  in  the  next  phase  of  the  study,  utilising  a  six-degrees-of-freedom  (6- 
DOF)  flight  path  simulation  computer  code,  called  FLIPSIM  (Flight  Path  Simulation  of  Missiles) 
and  assessment  of  systems  damage  and  functional  degradation. 

6  Discussion 


The  code  MISDAC  which  has  been  described  in  this  report,  is  a  derivative  of  TARVAC,  but  in¬ 
cludes  many  features  of  MIS  VAC.  Basically,  the  geometrical  target  descriptions  are  the  same, 
but  weapons  effects  options  have  been  reduced  (no  shaped  charge  effects,  for  example). 
However,  blast  damage  and  direct  hit  damage  from  intercepting  air-defence  missiles  have  been 
included,  with  the  general  assumption  that  blast  damage  does  not  prevail  the  overall  damage  as¬ 
sessment  of  missiles.  The  probability  of  kill  given  a  hit  has  been  assessed  by  two  independent 
probability  functions,  viz.  probability  of  damage  given  a  hit  and  probability  of  kill  given  damage. 
This  approach  simplifies  the  input  and  reduces  the  number  of  variables  compared  with  TARVAC 
while  still  retaining  some  “fuzzy”  approach  in  the  overall  assessment. 
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7  Conclusions 


A  missile  damage  assessment  code  MISDAC  has  been  derived  from  two  formerly  existing  target 
vulnerability  codes,  i.e.  the  "‘standard”  target  vulnerability  assessment  code  TARVAC  and  the 
derived  missile  vulnerability  assessment  code  MISVAC. 

The  final  kill  probability  assessment  still  has  to  be  implemented  using  advanced  six-degrees-of- 
freedom  simulations  to  establish  post-damage  missile  flight  trajectories,  which  determine  the 
overall  missile  kill  probability  given  engagement  by  an  air  defence  weapon.  The  new  code  as¬ 
sesses  the  damaging  effects  of  fragments,  projectiles  and  blast  utilising  probability  functions  and 
logical  expressions.  Fragmentation  damage  assessment  has  been  improved  by  implementing  ran¬ 
dom  variation  in  impact  shape  (number),  the  variance  being  dependent  on  shape  number  (sphere 
^  variance=0).  Blast  damage  assessment  has  been  simplified  using  blast-critical  lines 
representing  planes  and/or  components.  An  overall  blast  vulnerability  equation  employing  three 
“Westine”  constants  determines  the  relation  between  critical  distance  and  charge  weight  (R-W). 
Direct  hit  damage  is  also  incorporated,  which  in  general  is  assumed  catastrophic  for  the  target. 
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